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Abstract We use fluctuating hydrodynamics to evaluate the enhancement of ther- 
mally excited fluctuations in laminar fluid flow using plane Couette flow as a rep- 
resentative example. In a previous publication [J. Stat. Phys. 144 (2011) 774] we 
derived the energy amplification arising from thermally excited wall-normal fluctu- 
ations by solving a fluctuating Orr-Sommerfeld equation. In the present paper we 
derive the energy amplification arising from wall-normal vorticity fluctuation by 
solving a fluctuating Squire equation. The thermally excited wall-normal vorticity 
fluctuations turn out to yield the dominant contribution to the energy amplifica- 
tion. In addition, we show that thermally excited streaks, even in the absence of 
any externally imposed perturbations, are present in laminar fluid flow. 

Keywords Energy amplification ■ Fluctuating hydrodynamics • Laminar fluid 
flow ■ Orr-Sommerfeld equation ■ Plane Couette Flow ■ Squire equation ■ Vorticity 
fluctuations 



1 Introduction 

The presence of gradients, such as temperature gradients, concentration gradients, 
or velocity gradients, always causes non-equilibrium enhancements of thermal fluc- 
tuations that are spatially long ranged, even when the system is far away from any 
hydrodynamic instability pUJ. The present paper is part of a detailed study of the 
nature of thermally excited fluctuations in laminar fluid flow using plane Couette 
flow as a representative example. It has been verified that fluctuating hydrod> 
namics, originally developed for thermal fluctuations in equilibrium states 0, 
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can be extended to deal with thermal fluctuations in non-equilibrium states |2|. 
In fluctuating hydrodynamics the usual deterministic hydrodynamic equations are 
supplemented with random dissipative fluxes of thermal (natural) origin. In the 
case of laminar fluid flows one has to consider a random stress tensor to account 
for intrinsic thermal noise that always will be present. This noise will be ampli- 
fied by the presence of a velocity gradient. Energy amplification induced by the 
flow has attracted the attention of many investigators. In many of the studies 
reported in the literature the noise is not of thermal origin and does not obey a 
fluctuation-dissipation relation [5-9]. Fluctuating hydrodynamics provides a sys- 
tematic method for assessing the nature of spontaneous fluctuations in laminar flow 
induced by intrinsic noise. The application of fluctuating hydrodynamics to shear 
flows has been initiated by some previous investigators, but without considering 
confinement effects [10l4l4| . However, the long-ranged nature of the fluctuations 
is highly anisotropic and for certain directions of the wave vector the fluctua- 
tions encompass the entire fluid system, so that boundary effects will affect these 
fluctuations. 

Previously we have derived the appropriate fluctuating hydrodynamics equa- 
tions for laminar fluid flow. Specifically, we have shown how the thermally excited 
wall-normal velocity fluctuations can be described by a stochastic Orr-Sommerfeld 
equation [15] and the thermally excited vorticity fluctuations by a stochastic Squire 
equation [16j |. Accounting for realistic boundary conditions we obtained solutions 
of the stochastic Sommerfeld and Squire equations based on semi-quantitative 
Galerkin approximations [l5Hl7| . We now have derived more exact solutions of 
these stochastic equations in terms of an expansion of the eigenfunctions (hydro- 
dynamic modes) of the hydrodynamic operator. The more exact solution of the 
stochastic Orr-Sommerfeld equation for the wall-normal velocity fluctuations has 
been presented in a previous article in this series, to be referred to as paper I [18]. 
We found that the flow-induced enhancement of the wall-normal velocity fluctua- 
tions and the resulting energy amplification increases with the Reynolds number 
Re approximately as Re 2 . The actual enhancement of the velocity fluctuations 
strongly depends on the wave number. For large wave numbers (in the bulk of the 
fluid), this enhancement varies as the fourth power of the inverse wave number, in- 
dependent of any boundary conditions. For small wave numbers the enhancement 
vanishes as the square of the wave number due to the presence of boundaries. Our 
previous approximate solution based on a Galerkin approximation [l5| did repro- 
duce the correct dependence of the non-equilibrium enhancement of the fluctua- 
tions on the wave number but underestimated the magnitude of the enhancement 
at intermediate wave numbers [18( . 

The present paper is concerned with an analysis of the solution of the stochas- 
tic Squire equation for the wall-normal vorticity fluctuations. We shall proceed as 
follows. In Sect. [2] we recall the expressions for the stochastic Orr-Sommerfeld 
and Squire equations in terms of suitable dimensionless variables. The Squire 
equation for the wall-normal vorticity fluctuations includes a coupling with the 
Orr-Sommerfeld equation for the wall-normal velocity fluctuations. Hence, the so- 
lution of the stochastic Squire equation for the wall-normal vorticity fluctuations 
to be obtained in the present paper will depend on the solution of the stochas- 
tic Orr-Sommerfeld equation for the wall-normal velocity fluctuations obtained 
in our preceding paper (l8j ]. In Sect. [3] we describe the procedure for solving the 
stochastic Squire equation. For this purpose we expand the solution in terms of 
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the eigenfunctions of the linear hydrodynamic operator associated with the Squire 
equation. In Sect. 0] we derive the corresponding hydrodynamic modes and de- 
cay rates. In Sect. [5] we then obtain the expressions for both the equilibrium and 
nonequilibrium contributions to the intensity of the vorticity fluctuations, while 
in Sect. [6] we deduce the nonequilibrium energy amplification arising from these 
vorticity fluctuations. The nonequilibrium energy enhancement turns out to be 
proportional to the square of the Reynolds number. Specifically, we evaluate the 
nonequilibrium energy amplification associated with the vorticity fluctuations for 
fluctuations with wave vector in the spanwise direction which appears to be the 
most interesting case. We also find good agreement between the exact solution, 
obtained in this paper, and the semi-quantitative solution previously obtained in a 
Galerkin approximation [16j . We conclude with some general comments in Sect. [7] 

2 Fluctuating Hydrodynamics of Shear Flow. The Stochastic 
Orr-Sommerfeld and Squire Equations 

We consider a liquid with uniform temperature T under incompressible laminar 
flow with uniform density p between two horizontal plates separated by a distance 
2L. As in our previous publication [18] we adopt a coordinate system with the 
X-axis in the streamwise direction, the F-axis in the spanwise direction, and the 
Z-axis in the wall-normal direction [19| . Thus the mean flow velocity vn = {-yz, 0, 0} 
is in the X-direction with 7 representing a constant shear rate in the Z-direction. 
The bounding upper plate, at the position z = +L, moves in the positive X- 
direction with velocity jL, while the lower plate, at the position z = —L, moves 
with the same velocity in the opposite direction. This flow configuration is com- 
monly referred to as plane Couette flow. It is convenient to use a dimensionless 
position variable r, measured in terms of the length L, a dimensionless time t 
obtained by multiplying the actual time with the shear rate 7, a dimensionless 
fluid velocity v in terms of the product 7L, and a dimensionless stress tensor II 
in terms of pL 2 -y. 

We want to study velocity fluctuations around the stationary flow solution 
of the Navier-Stokes equation. Specifically, we are interested in fluctuations of 
thermal origin, i.e., fluctuations resulting from the intrinsically stochastic nature 
of molecular motions. Such fluctuations are always present and are unavoidably 
linked to any dissipative processes that are present in the system, Newton's law 
of viscosity in the present case. Fluctuating hydrodynamics provides a general 
and systematic framework for describing such thermal fluctuations, even when the 
system is in a non-equilibrium state 0, Q • The idea is that the linear phenomeno- 
logical laws representing dissipation in the system are to be supplemented with 
random dissipative fluxes (thermal noise), whose statistical properties are given 
by the fluctuation-dissipation theorem, see e.g. The goal is then to obtain 

the correlation functions of the fluctuating thermodynamic fields, velocity fluctu- 
ations <5v in our case, in terms of the statistical properties of the thermal noise, 
the stochastic stress tensor 511 in our case. Implementing this procedure we have 
shown in previous publications that the fluctuations 5v z of the wall- normal velocity 
component v z satisfy a stochastic Orr-Sommerfeld equation [lij 



d t (V 2 Sv z ) + z d x (V 2 Sv z ) - —X7 4 (Sv z ) = - {V x V x [V (<5JT)]} 2 , (1) 
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and the vorticity fluctuations Sw z = dxSv y — dySv x a stochastic Squire equation |l6| | 

d t {Sw z ) + z d x (5w z ) ~ d v 8v z - l-v 2 (Sw z ) = {v x [v {sn)}} z . (2) 

In these equations Re is the Reynolds number 

Re=^ (3) 

n 

with r\ being the shear viscosity of the fluid. By combining Eqs. ([TJ and ([2]) with 
the incompressibility assumption, V ■ <5v, one can obtain the three components of 
the fluctuating velocity field 6v. 

The difference between the stochastic Orr-Sommerfeld equation ([TJ) and the 
stochastic Squire equation ([2]) and their deterministic counterparts (lOL l2l| is the 
presence of noise terms on the right-hand side (RHS). These additive noise terms 
appear as derivatives of the random stress <5I7; their correlation functions can be 
deduced from the fluctuation dissipation theorem which in this case reads Q 

(M7«(r,t) ■ Sn kl (r', t')) = 2S (fofy + S u 5 jk ) <5(r - r') S(t - t'). (4) 

In this equation S is the dimensionless strength of the thermal noise [l8| : 

- k B T rj k B T 1 1 , . 

-y3 L 7 p 2 po L3^ L 2 Re - K°J 

where k B is Boltzmann's constant. We note that the actual correlation function for 
the fluctuating stress tensor only depends on the properties of the fluid, namely 
temperature, density, and viscosity; the shear rate and the Reynolds number only 
appears in (J5j> as a consequence of the manner in which the stress has been made 
dimensionless. 

Equations ([1]) and ([2]) form a pair of coupled stochastic differential equations 
which have to be solved for the velocity and vorticity fluctuations subject to ap- 
propriate no-slip boundary conditions: 

Sv z (r, t) = d z Sv z (r, t) = Sw z (r, t) = 0, at z = ±1. (6) 

In our earlier work we have accounted for these boundary conditions by solving 
Eqs. (fTJ) and ([2|) using a Galerkin method that allowed for analytical but approxi- 
mate expressions for the correlation functions of the wall-normal velocity and vor- 
ticity fluctuations [l5l[l6j . We have summarized these ap prox imate solutions of the 
two stochastic equations in a subsequently publication [17j . While the Galerkin- 
approximation technique enabled us to obtain relatively simple analytical results, 
the approximation is somewhat uncontrolled. Hence, we found it desirable to com- 
pare the approximate analytical solutions with exact numerical solutions that can 
be obtained through expansions in eigenfunctions of the hydrodynamic opera- 
tors. We have implemented this project for the solution of the Orr-Sommerfeld 
stochastic equation ([1]), yielding the correlation function for the wall-normal ve- 
locity fluctuations in paper I in this series jl8( . In the present paper we analyze 
the solution of the stochastic Squire equation ([2]) for the wall-normal vorticity 
fluctuations. There is an important difference between the Orr-Sommerfeld equa- 
tion ([TJ) and the Squire equation (J2j) . The solution of the Orr-Sommerfeld equation 
for the velocity fluctuations is independent of the solution of the Squire equation 
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for the vorticity fluctuations. On the other hand, the Squire equation ([2]) is cou- 
pled with the Orr-Sommerfeld equation though the presence of the term d y Sv z 
in Eq. Physically it means that the wall- normal velocity fluctuations involve 
only a coupling between the same (viscous) hydrodynamic mode (at different wave 
numbers) to which we have referred as " self-coupling" [l5| . However, the solution 
of the Squire equation for the wall-normal vorticity fluctuations is related to two 
mode-coupling mechanisms: a self-coupling between vorticity fluctuations and a 
cross-coupling between velocity and vorticity fluctuations. We shall recover our 
earlier observation [16] that the self-coupling mechanism in the Squire equation 
determines the intensity of the vorticity fluctuations in equilibrium and that the 
cross-coupling mechanism in the Squire equation determines the nonequilibrium 
intensity enhancement of the wall-normal vorticity fluctuations. 

To put our present work in the context of previous investigations of the effect of 
externally imposed stochastic forcing on shear flows we may mention the following. 
In their original papers, Farrell and Ionannou 0, 0] and Bamieh and Dahleh Q 
introduced a stochastic forcing term directly into the RHS of the Orr-Sommerfeld 
and Squire equations ([1) and ([2]) instead of the thermal forcing incorporated by 
us. More recently, Jovanovic and Bamieh Q introduced the forcing in the more 
basic Navier-Stokes equation, resulting in a more transparent interpretation of 
the random terms as actual forces. Moreover, Jovanovic and Bamieh Q allowed 
for some flexibility in the spatial spectrum of the noise, distinguishing between 
structured and unstructured noise. However, in spite of this flexibility, thermal 
noise was not included as a particular case considered by the authors [8J]. Indeed, 
from the flutuation-dissipation relation it follows that the spatial spectrum of 
the resulting stochastic forcing in the RHS of Eqs. (fTJ) and d2|) is not contained in 
the expressions of Jovanovic and Bamieh 8]. Hence, the results obtained here will 
differ from previous investigations. 

In terms of the input-output nomenclature that, borrowed from the dynamics 
and control literature, has become popular lately [22|; we take as input the 
thermal noise and our output (through the linearized Navier-Stokes equations) is 
the wall-normal vorticity autocorrelation function. 



3 Procedure for Solving the Stochastic Squire Equation 

Just as the procedure used for solving the Sommerfeld-Orr equation [is| . to solve 
the Squire equation ([2]) we apply a Fourier transform in time and in the horizontal 
XY-plane parallel to the walls: 

[icj + H] ■ <ku z (w,q||,z) = iq y 5v z (u, , z) + Sz{u,<l\\, z), (7) 

where uj is the frequency of the fluctuations 8w z and q|| = {q x , q y } the correspond- 
ing wave vector in the plane parallel to the walls. In Eq. ([7]), H represents a linear 
hydrodynamic operator: 

The RHS of Eq. ([8| contains two stochastic forcing terms. The first one, iq y Sv z , 
accounts for a stochastic forcing originating from the wall-normal velocity fluc- 
tuations, which can be represented in terms of the exact solution of the Orr- 
Sommerfeld equation with thermal forcing derived in paper I |l8(. The second 
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forcing term S z is given by the Fourier transform, in time and in the XY-plane, 
of the combination of derivatives of the random stress S TI in the RHS of Eq. ([2]) : 

S z = id z [q x SII Z y - q y 8Il zx ] - q x 8II X y + qySIIy X + q x q y [Sn xx - SII VV ]. (9) 

Since the hydrodynamic operator is linear we find the solution of Eq. (0) by ex- 
panding in a set of right eigenfunctions, i?jv(qii, z), or hydrodynamic modes as: 

oo 

5w*(w,q||,z) = ^2 Gjv(w,q||) Rn(m\\,z), (10) 
JV=o 

where the hydrodynamic modes are the solution of: 

h ■ r n {<\\\ , *) = rjv(q|| ) R N (% , *) (ii) 

with _T/v(q||) being the corresponding eigenvalue or decayrate. In Eq. (|lip R^(q^,z) 
must satisfy the boundary conditions i£jv(qii,±l) = 0. In Eqs. (|10p -(|ll | ) we antic- 
ipated the fact, to be discussed in more detail in Sect.[H that the right eigenvalue 
problem of the Squire operator, Eq. has indeed an infinite numerable set of 

solutions. Next, to evaluate the coefficients Gjv(w,qii) of the series expansion (|10[) 
we use the property that the complex conjugates i2^ r (qii , z) are the left eigenfunc- 
tions of the Squire operator ([8|), with corresponding eigenvalues Ijy(q||). Indeed, 
by using the boundary conditions ([6|) it can be readily shown that the adjoint of 
the operator H is simply its complex conjugate and the left eigenfunction is just 
the complex conjugate of the right eigenfunction. As a consequence, the biorthog- 
onality [23fl condition reads: 



I Rm(<1\\,z) Rn(<1\\,z) dz = B N (q\ 



) <W (12) 



Equation (|12[) determines the normalization of the eigenfunctions. Next, we evalu- 
ate the coefficients G^(lu,^) in the expression ()10p for 8uj z {ui, qn , z) by substitut- 
ing Eq. (110p into Eq. ((?])■ The resulting expression is then projected onto Rm{s\\i z ) 
and using the biortogonality condition (|12p . one readily obtains: 



Gjv(w,qii) = — -, u . " , — yr, (13) 

11 Biv(^,q||)H + -riv(q||)J 

where Fjv(w,qii) are the projections of the RHS of Eq. (JTJ) onto the R^(qpz) 
functions, namely: 



F N = J Rxiqpz) [iq v 5v z (u,qpz) + S x (u,q\\,z) 



dz. (14) 



For the evaluation of the vorticity fluctuations, we need the correlation functions 
(Fjy(jjj, q||) ■ F^f(w , qj| )). These, in turn, can be deduced from the correlations: 

(Sv* z (uj,qpz) ■ 5v z {J ,q\,z')) = (2tt) 3 5(u - J) 5{q\\ - q\) C zz {u, q\\, z, z), (15a) 
{S* z {u,q\\,z) ■ S z (u',q'pz')} = (2tt) 3 &{u-J) 5(q|| -qj|) 

x 25 q\{q\ + d z d z >)5(z- z), (15b) 

(S* z (u,q\\,z) ■ 5v z {J ,q\,z')) = (5vi(u,qpz) ■ S z (u',q\\,z')) = 0. (15c) 
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Equation (|15ap represents the solution of the stochastic Orr-Sommerfeld equation 
obtained in paper I [l8| (see Eq. (24) in [l8j]), where explicit expressions for the 
function C Z z{pj,c^\\,z, z ) in terms of the hydrodynamic modes and decay rates 
of the Orr-Sommerfeld operator have been presented. We shall not repeat those 
expression here, although they shall be used in some of the following calculations. 
Equation (|15c[l is a consequence of the fact that the random noise terms in the 
stochastic Orr-Sommerfeld and Squire equations are uncorrelated, see Eq. (18) in 
Ref. 0. Finally, Eq. (pE) . which is new in this paper, is readily obtained from 
the definition ([9|) of S z (a;,q||, z) and the fluctuation-dissipation theorem for the 
random stress tensor. The expression for the prefactor S is given by Eq. (|5|). 

With the help of Eqs. (|15|) . the correlations {Fpf(u, qy) ■ JjwC^'j q||)) between 
the coefficients defined in Eq. (|14p . can be expressed as: 

CFjy(w,q||) ■ F M (cj / ,q||)} = [ S^fay ) + 3$m ( w > Q|| )] ( 2n ) 3 S{u - u') %|| - qy), 

(16) 

with mode-coupling coefficients 

25 gjj I dz R* N (qpe) \q\ - df^ R M (q\\,z), (17a) 

1y J J dz dz ' r *n(<1\\,z) C zz (u,q.\\,z,z') R M (q\\,z'). (17b) 

The first expression, (|17a|l . is obtained from Eqs. (| 14p and (|15b|) if one performs 
integrations by parts to move the derivatives in (|15b|) from the delta function to 
the hydrodynamic modes, and uses the boundary conditions J2^r(qii,±l) = 0. In 
Eqs. (|16p and (|17[) we have introduced superscripts (E) and (NE) to distinguish be- 
tween the two contributions, anticipating the fact (to be discussed at length later) 
that the first set of mode-coupling coefficients will contribute only to the equi- 
librium equal-time vorticity fluctuations, while the second set of mode-coupling 
coefficients (due to the coupling of wall-normal velocity and vorticity fluctuations) 
contains the nonequilibrium amplification (enhancement) of those fluctuations. 

Equation (|17a,|) for "^^f (*3||) can be further transformed by integrating by 
parts while making use of the eigenvalue problem (|11|) . to obtain an expression 
more useful for future use: 

f' 1 

s nL^\\) = ^ Re 4 [rUi\\) + r M{q\\)] J d£ R*N(q\\,0 R M (qpt), (18) 

where we renamed the integration variable as £. 



*-NM (^qil) 



4 Hydrodynamic Modes and Decay Rates 

The hydrodynamic modes Rpf(q^,z) of the operator H are obtained by solving 
Eq. pip with the appropriate boundary conditions, i.e., Jijv(q||,±l) = 0. In view 
of the definition Q of H, the general solution of Eq. (|lip can be expressed as 
a linear combination of Airy functions. The boundary condition at z = 1 can 
be easily accommodated by a convenient selection of the coefficients. Then the 
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solution of Eq. satisfying the boundary condition at z = 1 can be written in 
the form $2$: 

R N {z) = Bi [{q x Re) 1/3 {a N - i)] Ai [{q x Re) 1/3 {a N - b)] 

-Ai[(^Re) 1/3 (ajv-i)] Bi [( 9a; Re) 1/3 ( aw -iz)] . (19) 

In principle, the parameter in fj 19[) can be any complex number. The decay 
rate of the hydrodynamic mode (|19[) is expressed in terms of this parameter apf 
as: 

ZV(q||)= fe ajv(q||) + JL. (20) 

For the hydrodynamic mode (|19p to satisfy the boundary condition at z = — 1 we 
need to impose the condition: 

= Bi [(^Re) 1/3 ( aA r - i)l Ai \(q x Re) 1/3 (a N + i)l 

-Ai [(^Re) 1/3 (a A r-i)] Bi [( g;r Re) 1/3 ( a7 v + i)] . (21) 

Because of the oscillatory character of the Airy functions, Eq. (|21[) has an infinite 
numerable set of complex roots a^. This fact has been anticipated in Eq. (|10p . and 
the index N has been used throughout to distinguish among the various modes. 
In principle, the decay rates depend on the parallel wave vector qn = {qx,<ly} 
and on the Reynolds number Re. However, as a consequence of the structure of 
Eq. (|2ip . the decay rates only depend on the magnitude q\i of the wave vector q|| 
in the plane parallel to the plates and on an effective Reynolds number 

Re = Re cosy, (22) 

where <p is the azimuthal angle of the wave vector qn , measured with respect to 
the stream-wise X-direction. Hence, g^Re = g||Re. This simplification is commonly 
referred to as Squire symmetry [Tjj]. As was discussed in detail in the preceding 
paper [18( , the eigenvalues and the eigenfunctions of the Orr-Sommerfeld hydro- 
dynamic operator obey the same Squire symmetry. However, in contrast to the 
decay rates of the Orr-Sommerfeld hydrodynamic operator, the decay rates ajv of 
the Squire hydrodynamic operator have an additional symmetry property, namely 

ajv(<7||,Re) = ajv(?||A,ReA) (23) 

for any real parameter value A. As a consequence, it is sufficient to determine the 
solutions of Eq. (|2ip for a single value of the effective Reynolds number, such as 
Re= 1. 

In Fig. [T] we show, as a function of the wavenumber qu , the decay rates of 
the Squire hydrodynamic operator H for Re = 1. Use of the explicit exact ex- 
pression (|19[) for the eigenfunctions (and, hence, Eq. (|21|l ) enables us for a much 
simpler computation of the data displayed in Fig. [T] when compared to a direct 
numerical integration of Eq. (|11[) performed by Gustavsson and Hultgren [25| . As 
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Fig. 1 Real (upper panel) and imaginary (lower panel) parts of the eight slowest decay rates 
as a function of gy, for effective Reynolds number Re = 1. 

expected, the results are the same: see, e.g., Fig. 1 of Ref. (25[. For small the de- 
cay rates are real numbers and a simple perturbative calculation allows to obtain 
the first terms in a series expansion in powers of 311 , namely: 

jV 2 7T 2 

a N{%)\ m=l = ^- + 0{ % ), for N= 1,2,3... (24) 

For larger qu and depending on the order N, the decay rates merge in pairs of 
complex conjugate numbers. For even larger — > 00 the real part of aj\[ decay to 

— 1/3 

zero as q,, , while the imaginary parts of each pair approach ±i. As explained 
above, this general landscape of decay rates is maintained for any Reynolds num- 
ber because of the scaling relation (|23|) . In particular, none of the decay rates 
becomes zero for any value of the wave vector or the Reynolds number. The same 
is true for the Orr-Sommerfeld equation and, as is well known, there is no linear 
hydrodynamic instability in plane Couette flow. 

One difference with the landscape of eigenvalues of the Orr-Sommerfeld oper- 
ator is that there is no transient merging of eigenvalues in a limited "window" of 
wave numbers, nor crossing of eigenvalues of different order. Hence, the curious 
behavior of the decay rates associated with the Orr-Sommerfeld hydrodynamic 
operator, as illustrated in Fig. 2 of paper I [l8| or in Fig. 2 of Gustavsson and 
Hultgren [25j , is absent in the case of the decay rates associated with the Squire 
hydrodynamic operator. The general nature of the eigenvalue map for the Squire 
operator is the same as depicted in Fig. [1] here, independent of the Reynolds num- 
ber. 

In view of the structure (|19p of the eigenfunctions and of the boundary condi- 
tions iljv(<Z|i,±l) = 0, the normalization constants (|12|) can be evaluated exactly 
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by using both the known integrals of products of the Airy functions [26( and the 
Wronskian of Airy functions. This procedure yields: 



B N (g\\) 



L _j Ai 



(g a Re) 1 /3( ajv _ i )' 



(feRe) 1 /3(a Ar + i ) 



(25) 

and a similar expression in terms of the Bi functions, see Eq. (|21[) . Again, as 
a consequence of the Squire symmetry, the normalization constant Bj^(q^) only 
depends on the product g x Re = g||Re. 



5 Nonequilibrium Vorticity Fluctuations 

Starting from the general theory presented in Sect. [3] and using the expressions 
for the hydrodynamic modes and decay rates derived in Sect. [3] we can now evalu- 
ate the autocorrelation function (Sw z (uj,q^, z) ■ 8w z (lj', q|| , z')) of the wall-normal 
vorticity fluctuations. By combining Eqs. (|10p . ()13[) . and (1161) . one readily obtains: 

(5w* z (u>,q\\,z) ■ Sw z (uj',q'pz')} = W zz (uj, q\\ , z, z) (2tt) 3 5(ui - u) 5(q\\ - q'|) (26) 

with 

°° ^(E) , -(NE), \ 

W zz (u>, n ,z,z>)= £ — -piM_t "ATM} ) RMz)Rm{z >). ( 27) 

To simplify the notation we have suppressed on the RHS of Eq. (|27p the explicit de- 
pendence of the quantities on the wave vector qii . In principle, Eq. (|27[) enables us 
to investigate not only the intensity of the vorticity fluctuations, but also the time- 
dependent correlation function characterizing the dynamics these fluctuations. 
However, just as in the case of the velocity fluctuations derived in paper I |l8|, we 
consider here only the intensity of the nonequilibrium vorticity fluctuations, which 
is given by the equal-time correlation function (6u) z (q\\,z,t) ■ Sw z (q'^, z' ,t)}. This 
equal-time correlation function is obtained by applying a double inverse Fourier 
transform, in frequencies uj and uj' , to Eq. (126[) so that 

(Scj* z {qpz,t) ■ 6w z {q\\,z',t)} = W zz (q\\ , z, z) (2tt) 2 5(q|| - qj|), (28) 

with 

1 f°° 

Wzz{m,z,z') = — du;W zz (uj,m,z,z'). (29) 

^ J-oo 

Substituting Eq. (|27|) into Eq. (|29|). and performing the integration over the fre- 
quency uj of the fluctuations, one readily obtains an expression for the amplitude 
of the equal-time correlation function. In view of the structure of Eq. (|27[) , we con- 
clude that the resulting expression will contain two additive contributions, namely 
an equilibrium contribution (E) and a nonequilibrium contribution (NE): 

W zz (q,| , z, z') = Wif (q,| ,*,*') + W4 NE) (q,| ,Z,z'). (30) 
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5.1 Equilibrium contribution to the intensity of the fluctuations 

Since the mode-coupling coefficients £\tm ^° no ^ depend on the frequency w, the 
integration in Eq. (|29[) can be readily performed for the equilibrium contribution. 
Combining the result with the expression (|18p for the "equilibrium" mode-coupling 
coefficients we arrive at: 

W £\ % ,z,z')=~SKz4 f; f dt R »® l Mi0 R%(z) R m (z'), 

Af,JU=0 1 N M 

= SRe q\ £ f dC^Mil R* N (z) 6(t /), (31) 
N=Q J-i v N 

= SRe qu 8{z — z), 

where we have made use of the expansion of the delta function in terms of the 
hydrodynamic modes Rpf{<\\\, z), 

oo 

s (z - *) = E wn Rn ^ ' c) Rn ^ ' z) ' (32) 

n=o N[ - q V 

which is obtained by using the orthogonality condition (|12[). It is obvious that 5(£ — 
z), as a function of z, satisfies the relevant boundary conditions. Expression (|32[l is 
valid for any value of the wave number qu, or the effective Reynolds number Re. The 

important result is that the expression (|3ip for wjf^ (qii, z, z') indeed reproduces 
the intensity of the equal-time autocorrelation function of the wall-normal vorticity 
fluctuations for a fluid in equilibrium. Note that the prefactor 5Re appears in 
Eq. (|3ip as a consequence of the adoption of dimensionless variables; when one 
reverts to physical variables the resulting prefactor is indeed independent of the 
shear rate. 

Equation (|31[) justifies the use of superscripts (E) and (NE) in Eq. (|16[). In ad- 
dition, it shows that the nonequilibrium contribution to the intensity of vorticity 
fluctuations arises only from the coupling with the wall-normal velocity fluctua- 
tions in the stochastic Squire equation ([2]), and not from the self-coupling also 
present in the stochastic Squire equation ([2]). This result was already found pre- 
viously [3 on the basis of a Galerkin approximation. We now see that this is 
exact, and not just a consequence of the simplicity of the approximation used pre- 
viously However, the self-coupling in the Squire contribution does contribute 
to the time-dependent nonequilibrium correlation function which is not considered 
here. 

The separation (|30p of the effects of thermal noise into an equilibrium and 
a nonequilibrium contribution is equivalent to what was found in the study of 
stochastic forcing, for instance, Eq. (17) of Ref. [Tj or Eq. (4.3) of Jovanovic and 
Bamieh [8| . The terms proportional to the cube of the Reynolds number in those 
expressions @, H| vanish when there is no flow in the system, but physically, the 
terms proportional to Re exists even when there is no flow. Indeed, as is the case 
in our Eq. (|3ip . one power of Re in Refs. @, 01 is due to dimensionality reasons. 
A forcing (stochastic or not) must have units of force and the dimensionless time 
is in unit of the (inverse) shear rate. An important feature of our calculation is 
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that we recover the well-known expression (|31[) for the vorticity fluctuations in 
equilibrium [27L [28| . 



5.2 Nonequilibrium contribution to the intensity of the fluctuations 

Upon substitution of the (NE) part of Eq. (p7|) into Eq. ([29)1 . one obtains the 
nonequilibrium contribution Wz^ (qy , z, z') to the vorticity fluctuations. By mak- 
ing further use of Eq. (|17p for the mode-coupling coefficients, we arrive at the 
explicit expression: 



N,M=0 N 



1: 



In principle one needs to substitute the autocorrelation function C Z z(w, £, f') of the 
wall-normal velocity fluctuations from Paper I [18] into Eq. (|33p . and perform the 
various integrations and summations to obtain the function of interest. This pro- 
cedure yields rather complicated expressions and only marginal analytical progress 
can be made, at the expense of very large and cumbersome expressions. There- 
fore, unlike our previously solution obtained on the basis of a Galerkin approx- 
imation |l6j . only a numerical computation of Wz^ E \z, z') is generally possible, 
and even this numerical procedure turns out to be rather long and difficult. 

There is one important issue that should be mentioned. Because of the pres- 
ence of the term as a prefactor in Eq. (|33p . it turns out that the intensity of 
the nonequilibrium vorticity fluctuations has a maximum in the spanwise direc- 
tion (q x = 0), while it is zero in the streamwise direction. This is opposite to 
the wave-vector dependence of the wall-normal velocity fluctuations discussed in 
Paper I [l8j . In that case the intensity of the fluctuations has a maximum in the 
streamwise direction (q y = 0) and is zero in the spanwise direction. Moreover, for 
the same Reynolds number the intensity of the vorticity fluctuations is substan- 
tially larger than the intensity of the wall-normal velocity fluctuations. Hence, an 
important conclusion that can be derived from Eq. (|33|) is that the most impor- 
tant effect of the flow on thermal fluctuations is the enhancement of wall-normal 
vorticity fluctuations with wave vector in the spanwise direction; or, equivalently, 
fluctuations that are constant in the streamwise direction. The same conclusion 
was obtained from our previous approximate Galerkin solution [THl H(jj . see in 
particular Fig. 6 in Ref. [f?|. A similar conclusion is obtained both from direct nu- 
merical simulations of the full Navier-Stokes equations or from analytical studies 
of transient growth (amplification) of perturbations, see, for instance refs. (2^.[29j 

Not only are the vorticity fluctuations in the spanwise direction the most domi- 
nant and interesting. In addition, further analytical progress is also possible for this 
case. This is the reason why previous investigators, who considered externally im- 
posed forcing, have focused on vorticity response in the spanwise direction 0,0, EJ. 
Indeed, when q x = the (Fourier transformed) Orr-Sommerfeld and Squire equa- 
tions simplify notably. The eigenfunctions of the Squire operator can be simply 
expressed in terms of trigonometric functions, while the corresponding eigenvalues 
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can be obtained analytically. The eigenfunctions of the Orr-Sommerfeld opera- 
tor in this limit can be written as combinations of trigonometric and hyperbolic 
functions, while the eigenvalues can be obtained numerically by solving relatively 
simple algebraic equations, as first discussed by Dolph and Lewis [13] ■ One impor- 
tant property is that in both cases the eigenfunctions possess a well-defined parity, 
and can be naturally classified into odd eigenfunctions and even eigenfunctions. 
We note that this is not true in the general ^ 0. Hence, in the remainder 

of this paper we shall focus on the vorticity fluctuations in the spanwise direction, 
and their effect on the nonequilibrium energy amplification induced by the fluid 
flow. 



6 Nonequilibrium Energy Amplification 

We recall that in the derivation of the stochastic Orr-Sommerfeld and Squire equa- 
tions we have used the incompressibility assumption V • <5v = [l5l llq ] . Hence, 
only two components of the velocity fluctuate independently. Most investigators 
on the subject have been interested in the so-called kinetic-energy amplification, 
that can be obtained from the sum of the equal-time autocorrelation functions 
(5v* (qy , z, t) ■ dVi{<iH,z ,t)}. And indeed, as elucidated in more detail in our previ- 
ous analysis on the basis of a Galerkin approximation [16j, the spatial spectrum 
of the kinetic-energy amplification is proportional to the vertical average 



i j j dz dz' y^(foj(q||,z,t) •5v i {q[\\,z',t)) 



(34) 



dz dz I —{8w*(c^,z,t) ■ 8w z {q[n,z ,t)) + (5vt(^,z,t) ■ 5v z {q[»,z ,t)) 
l J-i 



i r 1 f\ . > i i 

2 



where on the RHS of this equation we only need to consider the wall-normal 
velocity and vorticity fluctuations because of the divergence-free condition V ■ <5v = 
0. In addition, the fact that there is no cross-correlation between Sw z and 8v z has 
also been employed. 

The contribution in Eq. (|34[) arising from the wall-normal velocity fluctuations 
has been studied extensively in Paper I [l8j|- We focus here on the contribution 
from the vorticity fluctuations which, in view of Eqs. (|28|) - (|29l) . will be proportional 
to the quantity 

W zz {^) = ^ J J dzdz' W zz (d {{ ,z,z'), (35) 

which is a function of the horizontal wave vector qy of the fluctuations that we 
shall investigate in this section. 

First of all, one notes that because of the structure of the mode-coupling coef- 
ficients, the spectrum W / Z z(q||) can be expressed as the sum of an equilibrium and 
a nonequilibrium contribution which we prefer to write in the form 

Wzz (q,| ) = W4 E) (?|| ) [l + AW^ E) (q,| )] , (36) 

where z^wj^^qn) represents the nonequilibrium energy enhancement. The equi- 
librium contribution W-uP {q\\ ) in Eq. (134[) is obtained by substituting Eq. (|31|) into 
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Eq. ([35]) 

Vl4 E) (9||) = SRegf, (37) 

which is the same as in the absence of any flow. 

The nonequilibrium enhancement AW^^qii) in Eq. (|36[) , resulting from the 
wall-normal vorticity fluctuations, can be obtained by substituting Eq. (|33f) into 
Eq. (|35[) . As already mentioned before, in general, this procedure yields a compli- 
cated expression that can only be evaluated numerically. 



6.1 Enhancement of streamwise-constant fluctuations (i.e. 
the spanwise direction) 



with wave vector in 



A particular simple case is that of fluctuations constant in the streamwise direc- 
tion with q in the spanwise direction, i.e., for which q x = and q v = q«, = q. 
As mentioned above, under this condition the working equations simplify greatly, 
and a more compact analytical expression can be obtained for the nonequilibrium 
enhancement. Because of its obvious simplicity,this particular case has been ana- 
lyzed in detail by some previous investigators [M S but for externally imposed 
stochastic forcing, not thermal noise. Hence, their results differ from the ones ob- 
tained here. For this reason, we present our explicit results for the enhancement of 
vorticity fluctuations induced by thermal noise with the (horizontal) wave vector 
in the spanwise direction, namely: 



1 9^=0 



9 — tanh q tanhq 9 tanh q 



2q 



cosh q 



(2^ + 3^) tanh x /a 2 N + 2q 2 



2q 2 

q 2 Ja 2 +2q 2 



+ Re 2 J2 



cosh 2 \ a 2 N + 2q 2 



N=0 



(a 2 N + q 2 ) 2 (a 2 N +2q 2 ) 2 



q tanh q 
< + 2q 2 



tanh » / a 2 N + 2q 2 



(38) 



with 



a N 



A2N + 1), 



(39) 



so that (a% + q 2 )/Ke are the eigenvalues of the Squire operator in the spanwise 
direction with corresponding eigenfunctions of even parity. As already mentioned, 
in this particular case (spanwise direction) the hydrodynamic modes (eigenfunc- 
tions) have a well-defined parity. Because of the ^-integrations in Eq. (|17b|) for the 
mode-coupling coefficients, modes with different parity do not couple. In addition, 
because of the integrations in Eq. (135[) . only the even eigenfunctions or modes 
do finally contribute to the nonequilibrium enhancement. To obtain Eq. (|38[) we 
have closely followed a procedure used by Bamieh and Dahleh [7]; in particular 
we used the auxiliary function g(z) introduced in their Lemma 4- However, we 
obtain a different result because the thermal noise, considered here, has a special 
spatial spectrum given by the fluctuation-dissipation theorem that is differ- 
ent from the spectrum of the externally imposed stochastic forcing considered by 
Bamieh and Dahleh Q . One power of the Reynolds number appears in the energy 
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q 

Fig. 2 Nonequilibrium enhancement AW^ E (q) of wall- normal vorticity fluctuations in the 
spanwise direction (q x = 0) from Eq. (138 6 - The two straight lines represent the asymptotic 
behaviors at small and large wave numbers, Eqs. (1406 and 14U . respectively. Data are for 
Re = 500. 



amplification because of dimensionality reasons. Physically, the flow-induced am- 
plification of thermal noise is proportional to the square of the Reynolds number 
(shear rate), not to the cube of the Reynolds number as stated by Bamieh and 
Dahleh [7] or Jovanovic and Bamieh Q . 

In Fig. [2] we show a plot of the nonequilibrium enhancement AW^(q) of 
the wall-normal vorticity fluctuations with wave vector in the spanwise direction 
{qx = 0, or streamwise constant) as given by Eq. (|38|) . The data in Fig. [2] are for 
Re = 500, but we note that the Reynolds number appears in Eq. t|38|) only as 
as a prefactor, so that the ratio zl W^ E (q)/Re does not depend on the Reynolds 
number. We conclude from Fig. [2] that the spanwise energy amplification can be 
very well visualized as a simple crossover between the two asymptotic behaviors, 
at large and small q, that can be easily obtained from Eq. (|38|) . namely 

&W ZZ (q) > Re q 



Q.r 



34 256 2 + cosh(2JV + l)7r 



ZOO \ - 

315 ~ ^ Q (2iV+ l) 7 sinh(2iV + 1)tt 
14 x 10~ 3 ReV, (40) 



Re 2 



AW™(q)\ (41) 
\q x =0 £q 

We recover in Eq. (|41[) the q -4 behavior that is typical of nonequilibrium fluctu- 
ations [2] at large wave numbers corresponding to wavelengths smaller than the 
spacing between the walls. For such wavelengths the fluctuations are not affected 
by the boundary conditions and we recover, as an asymptotic limit for large q, 
results obtained by previous investigators for nonequilibrium fluctuations in fluids 
under shear [10h14| . On the other hand, the vanishing of the intensity of fluctua- 
tions with small q, as implied by Eq. (140[) . is to be expected in the sense that the 
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walls (boundary conditions) effectively suppress fluctuations of very long wave- 
length, comparable with the separation distance between walls. From Fig. [2] we 
note that the flow at Re = 500 causes an enhancement of the thermal energy 
about 1000 times larger that the thermal energy that would be expected from a 
local equilibrium assumption. Such a profound enhancement of the fluctuations is 
a general phenomenon in fluids in nonequilibrium states Q ■ 

Figure [2] shows that the main effect of the flow on the thermal fluctuations is to 
select and maximally amplify the wall-normal vorticity fluctuations with a partic- 
ular value of the wave vector q m . As already discussed, the wave vector maximally 
enhanced is in the spanwise direction, and from Eq. (|38|) we find numerically its 
magnitude to be q m ~ 1.4103, which is the location of the maximum in Fig. [2j 
Therefore, in real space, the thermal noise amplified by the flow will manifest 
itself mainly as a set of vortices of size ~ (2tt/1A1)L ~ 4.5L distributed in the 
spanwise direction and that are constant (extremely elongated) in the streamwise 
directiorQ. The intensity of these fluctuating vortices will be proportional to k^T 
and to the square of the Reynolds number, see Eq. (J38J . But we note that close 
to the center of the layer, where the base flow velocity is zero, the velocity fluctu- 
ations might be of the same order as the mean (base) velocity. These fluctuating 
vortices deform the base flow, developing a series of streaks, i.e., narrow regions 
where the streamwise velocity is larger or smaller than the average, as indicated in 
Fig. [3] where we have tried to illustrate schematically this physical situation. These 
streaks will be typically separated by a distance of about 4.5 times the half gap 
between the walls, as also shown in Fig. [3] This optimal spanwise wave number is 
independent of the Reynolds number. 

At this point it is interesting to note that several authors have identified the 
appearance of a set of fluctuating streaks in sheared flows as a precursor of the 
instability (3lT l32l | . More quantitatively, extensive numerical simulations of the 
unstabilization of plane Couette flow |22i . [29, 33, 34J have revealed large-scale co- 
herent streaks with typical spanwise wavelength of » 3L — AL. For instance, most 
recently Gayme et al. [29j performed an extensive analysis of direct numerical sim- 
ulation data by Tsukahara et al. [35J], obtaining an optimal spanwise wavelength 
of 1.8 times the distance between the plates, or 3.6 times the half distance, to 
be compared with our result of 4.5L. Of course, the scope of our present work is 
restricted by the use of linear equations, so that its relevance to shear flow in- 
stability has to be considered, at most, as tentative. It is well-known and widely 
accepted that a complete understanding of unstabilization of shear flows requires 
a fully nonlinear theory [3lT I32I]. In any case, the identification of the modes that 
are maximally enhanced in a linear theory may provide useful insights for develop- 
ing simplified nonlinear theories, like the single mode model recently discussed by 
Gayme et al. |29l |. Furthermore, it is intriguing to know that a linear theory pre- 
dicts that thermal (natural) noise develops into "streaks", of intensity oc fexRe 2 , 
separated by a distance about 4.5L, in agreement with numerical simulations of 
the nonlinear problem [22I ]. 

To conclude this section, we mention that the algebraic oc g~ 4 / 3 dependence 
found for the enhancement of bulk fluctuations vanishes in the spanwise direc- 



1 This is by approximating the complicated spatial spectrum of the thermal velocity- 
fluctuations, as shown for instance in Fig. 6 of Ref. Il7ll . by just two delta functions located at 
the symmetric maxima (of the fluctuating wall-normal vorticity) at q x = 0, q y = ±1.4103. 
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Fig. 3 When we add to the mean Couette flow a set of wall-normal vortices with spanwisc 
modulation of about 4.5L and very elongated in the streamwise direction, narrow regions where 
the streamwise velocity is larger or smaller than the average develop, customarily referred to 
as streaks. For clarity, we draw these vortices only at a reference height, although in general 
they will present some vertical (wall-normal) extension. 



tion [l6(. That is the reason why the "shoulder" at intermediate q shown in 
the bottom panel of Fig. 2 of Ref. does not appear in our current results. 
Bulk nonequilibrium fluctuations in sheared fluids have been first investigated by 
Tremblay et al. who described the q~ A behavior. Subsequently, Dufty and 
Lutsko (ill Il3j included the algebraic wave-number dependence that, in general, 
appears at shorter wavelengths. Although these earlier papers 

0, El, El refer 

specifically to wall-normal velocity fluctuations, we have found elsewhere [16j that 
these features are also present for vorticity fluctuations in bulk. 



6.2 Comparison with Galerkin approximation 

In our previous publications jl6l Il7( the same problem considered here was in- 
vestigated on the basis of a simple Galerkin approximation. In these previous 
publications we did not present explicit analytical expressions for the intensity 
of the vorticity fluctuations (on the basis of a Galerkin approximation), but only 
displayed the results graphically. However, as in the case of the exact solution, 
also in case of the Galerkin approximation the solution simplifies greatly for vor- 
ticity fluctuations in the spanwise direction. Specifically, for the fluctuations in 
the spanwise direction the expression in terms of the Galerkin approximation in 
Refs. [H [13 reduces to 



W zz {q) = S | Reg 2 



1 + 



27ReV 



7(2g 2 -r-5)(4g 4 + 23g 2 + 78) 



(42) 
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0.1 1 10 
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Fig. 4 Comparison of the exact nonequilibrium enhancement AW^ E (q) of wall- normal vortic- 
ity fluctuations in the spanwise direction (q x = 0) from Eq. 1138 II (solid curve), with a Galcrkin 
approximation developed in previous publications 16, 17], sec Eq. i l-l-'li (dashed curve). Data 
are for Re = 500. The agreement is good. 

Equation (|42[) contains an equilibrium contribution that is about 17% percent 
lower than the exact result, Eq. (|37p . This is a shortcoming of the Galerkin- 
approximation scheme, as discussed elsewhere [l^]- But our purpose here is to 
compare the nonequilibrium enhancement, given by the second term inside the 
square brackets in Eq. (|42[l . with the exact result given by Eq. (|38|l . For this pur- 
pose we show in Fig [3] the exact AW^ E (q) in the spanwise direction from Eq. (|38p 
as a solid curve, together with the Galerkin approximation (|42[) of this quantity 
as a dashed curve. The agreement is rather good qualitatively. Quantitatively, 
the Galerkin approximation (|42[) for large q underpredicts the exact asymptotic 
limit (|41[) by 4% (indistinguishable on the scale of Fig. 2]) , whereas in the small 
q limit Galerkin approximation overpredicts the exact asymptotic limit (|42[) by 
about 60%. The prediction, on the basis of the Galerkin approximation (142[) . of 
the wave number of maximum enhancement q m is quite good, differing by about 
1% from the true value. 

Although we have here compared explicitly the exact solution with the Galerkin 
approximation only in the spanwise direction, from the good results obtained we 
may infer that the Galerkin approximation developed in Ref. will also yield 
a good representation of the nonequilibrium enhancement of the wall-normal vor- 
ticity fluctuations in any direction. This expectation has been confirmed by some 
preliminary calculations for arbitrary qn on the basis of Eq. (I33|) . 

7 Concluding remarks 

We have shown that even in the absence of any external perturbations, already 
the thermally excited fluctuations, which are always present, cause a substantial 
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energy amplification in laminar fluid flow with the main contribution to the en- 
ergy amplification arising from wall-normal vorticity fluctuations as a result of a 
coupling of these wall-normal vorticity fluctuations with the wall-normal velocity 
fluctuations. On the other hand, in computational fluid dynamics, to destabilize 
shear flows and to investigate the transition to turbulence, some externally imposed 
random initial conditions are customarily introduced [34| • We have demonstrated 
that such externally imposed initial conditions are not needed physically. Indeed, 
thermal noise causes already the base flow to spontaneously develop streaks, which 
are currently expected to be the first step in the transition to turbulence, even in 
a linear approximation. We have evaluated the typical spanwise distance between 
these thermally excited streaks, finding excellent agreement with what is observed 
in simulations. Of course, the problem of the nonlinear evolution of these streaks 
needs to be further investigated. For this purpose direct numerical simulations of 
fluctuating hydrodynamics may be of interest. 

In the last few decades there has been considerable interest in the fluid dynam- 
ics community in these and similar problems, like optimal disturbances, transient 
growth of perturbations, and energy amplification. At the same time, and some- 
what independently, in the statistical physics community there has been an interest 
in the nonequilibrium enhancement of thermal fluctuations. The purpose of this 
paper, as well of the previous paper [l8], has been to make a connection between 
the approaches used in statistical physics and fluid dynamics. 
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